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In extending the conventional dynamic models, we consider 
a simple model to account for the environment fluctuations of 
particle atoms in a protein system and derive the elastic inco- 
herent structure factor (EISF) and the incoherent scattering 
correlation function C(Q,t) for both the jump dynamics be- 
tween sites with fluctuating site interspacing and for the diffu- 
sion inside a fluctuating sphere. We find that the EISF of the 
system (or the normalized elastic intensity) is equal to that 
in the absence of fluctuations averaged over the distribution 
of site interspacing or sphere radius a. The scattering cor- 
relation function is C(Q,t) = J]„( e ~ A ' l( ' l)t > VK*)> where the 
average is taken over the Q-dependent effective distribution 
of relaxation rates A n (a) and ^(i) is the correlation function 
of the length a. When ^(j(t) = I, the relaxation of C(Q,t) is 
exponential for the jump dynamics between sites (since \ n (a) 
is independent of a) while it is nonexponential for diffusion 
inside a sphere. 

PACS numbers: 02.50.Ey ; 66.10.Cb ; 87.15.Hc 

By analyzing the incoherent scattering function, 
5(Q, uj), where Q is the scattering wavevector and Tiu> the 
energy transfer, techniques of quasielastic neutron scat- 
tering from hydrogen atoms (the main neutron scatterers 
in a typical protein) allow to study motions of particles 
(atoms, molecules, chemical species, ...) in biological sys- 
tems [GJ. For this purpose, there are many physical situ- 
ations of interests like molecules at surfaces, in micellar 
systems or in vesicles and structural cages, in which the 
system is modeled by jump of particles among sites or 
by diffusion inside a confining geometry ]^-||. For in- 
stance, when the problem can be described by the jump 
dynamics of a particle between two non-equivalent sites 
separated by a distance a, the elastic part oiS(Q,ui) (i.e., 
the elastic incoherent structure factor (EISF)), is given 
by 



A (Qa) = l-2p(l-p) [l-jo{Qa) 



(1) 



where ji(z) is the spherical Bessel function of the first 
kind of order / and p is the probability of finding the par- 
ticle in one of the sites. The incoherent scattering corre- 
lation function C(Q, t) (i.e., the inverse Fourier transform 
of the quasielastic part of S(Q, u>)) is a single exponential 
independent of Q given by C(t) = e~ r * where T (related 
to mean residence times of the particle in either site) is 
the relaxation rate of the probability of finding the par- 
ticle in either site. This jump model has been used to 



study the internal molecular motions and to interpret 
the dynamical transition in proteins 

Likewise, when the problem can be described by the 
isotropic diffusion inside a sphere with impermeable sur- 
face g , one finds that for Q small, the quasielastic part of 
S(Q, uj) is well described by a single Lorentzian with the 
linewidth T = 4.33L>/a 2 , where a is the sphere radius 
and D the diffusion constant of the particle. Equiva- 
lently, C(Q, t) ~ e~ r °* where T^ 1 is the typical time for 
a particle to diffuse over the entire sphere. For times of 
order or larger than Tq 1 (i.e., in the u> -C To limit) the 
EISF is given by 



Ao{Qa) = \F(Qc 



3ji(Qa) 



Qa 



(2) 



Such an analysis is used by several authors ]4|,|3||9| to 
study, for instance, the internal dynamics, structure and 
dynamics of surface molecules in proteins. In these stud- 
ies, the particles (mainly hydrogen atoms) are assumed 
to diffuse within permanent spherical cages. 

On the other hand, it is well known that proteins are 
fluctuating systems which undergo configurational fluc- 
tuations in their structure. In different conformational 
substates, a protein may have the same coarse structure 
but differs in local configurations leading to the fluctua- 
tions of the local environment of each protein atom. As a 
result, a distance between sites or a structural cage in a 
protein fluctuate in length or size (and shape) and have a 
finite lifetime due to local structural relaxation just like 
in mode-coupling picture of liquids ||To|| . In this respect, 
neglecting the fluctuations of shapes for simplicity, the 
dynamics of each hydrogen atom in a typical protein can 
be described, in the first approximation, as a jump dy- 
namics between sites separated by fluctuating distances 
or by diffusion inside a fluctuating sphere. 

In this paper, we present a simple model analysis to ac- 
count for fluctuations of the local environment. We focus 
on the derivation of the EISF and C(Q,t) in situations 
in which the site interspacing a (for the jump between 
sites) or the sphere radius a (for the diffusion inside a 
sphere) is allowed to fluctuate in the course of time. To 
model simply the fluctuations of a, we consider the fol- 
lowing dynamics for the system formed of a particle atom 
jumping between two sites (or diffusing inside a sphere 
with impermeable surface): the particle keeps jumping 
between sites separated by a(t) (or diffusing within the 
spherical cage of radius a(t)) until it suffer a configura- 



1 



tional collision of zero duration that equilibrates both the 
site interspacing (or the sphere radius) and the particle 
position. That is to say that after a collision the update 
particle position and the site interspacing (or the sphere 
radius) are chosen according to the normalized equilib- 
rium distribution, 



-/3[V(a)+U(r,a)] 



(3) 



where r = (r, O), r is the particle position, f2 stands for 
polar and azimuthal coordinates, U(r, a) is the potential 
energy for the particle position for a fixed a while V(a) 
is the potential of free energy associated to the length 
(spherical cage of radius) a, and /3 -1 = keT is the ther- 
mal energy. The reduced equilibrium distribution for the 
length a is defined as, p cq (a) — ' dQ, J^°r 2 P eq (r, a) dr. 
The waiting time between successive configurational col- 
lisions is a random variable with distribution <j){t) related 
to the stationary correlation function of site interspacing 
(or the sphere radius) by 



(a(t)a(O)) 



(-r)dr 



(4) 



Assuming that the environment fluctuations act on the 
particle dynamics in renormalizing just the residence 
times (for jump dynamics) or the diffusion constant (for 
diffusion dynamics), thus the Green's function, i.e. the 
probability density of finding the particle at the position 
r with site interspacing a (inside a sphere of radius a) at 
time t given that it was initially at i"o with site interspac- 
ing ao (inside a sphere of radius oq), is given by 



G(r, a,t\r ,a ) 



G (r,t\r ;a)S(a- a )ip(t) 
+ P cq (r,a) [l-if>(t)] , 



(5) 



where Go(r, t\rg; a) is the Green's function for a parti- 
cle with a fixed site interspacing a (or diffusing inside a 
sphere of a fixed radius a) . 

The incoherent intermediate scattering function, 
/(Q,i), (which is the inverse Fourier transform of 
S(Q,w))isI(Q,t) = 

/ da J dr J da J dr e lQ r G(r, a, t\r , a ) e" lQ r ° P oq (r , a ). 
Using into this relation the expression of the Green's 
function in Eq.(|J), one can show that for an isotropic 
problem the function J(Q,i) splits into two parts as, 

I(Q,t)=A(Q) + [l-A(Q)]C(Q,t), (6) 

where the elastic part, i.e. the EISF of the system, is 



A(Q) = 



poo pA-k poo 

/ da dVL r 2 e lQ r P cq {r 7 a)dr 
Jo Jo Jo 



(7) 



and the incoherent scattering correlation function C(Q, t) 
containing all information about the relaxation dynam- 
ics, is 



C(Q,t) = ^g„(Qa)e- A "( Q )y m 

°° p OO 

V / dag n [f3V(a),Qa]e- 
„_i Jo 



An(a)t 



m, (8) 



where the average of any function f(a) is defined as 
(/(«)> = fo^Pcqia) f(a) da. The g n {Qa) is the Q- 
dependent bare distribution of relaxation rates A„(a) in 
the absence of fluctuations while g n [(3V{a), Qa] is the Q- 
dependent effective distribution of relaxation rates which 
accounts for distribution of a. To be specific, we assume 
that the effective potential of free energy associated to 
the length a is given by 



V(a) = 



aa 

oo 



2k B T ln(a) 



a> R. 
a < R . 



(9) 



where e = oR 2 , 7? (of order of the Van der Waals contact 
length) is the minimum value of a, a is the force con- 
stant and the repulsive logarithmic term represents the 
entropic contribution which accounts for the increase in 
configuration space as a increases. As a result of the bal- 
ance between the attractive harmonic and repulsive en- 
tropic terms, the length a has an equilibrium value at fi- 
nite temperature T given by a cq = max [i?, (k^T/cr) 1 / 2 ] . 
We note in passing that such a harmonic potential in 
Eq. (^|) , without the entropic term, has been used to study 
the hydrophobic effect in the volume fluctuations of glob- 
ular proteins [jf]]] . 

In what follows we derive in detail the EISF and the 
C(Q, t) in Eqs.(|7]) and (||), respectively, for the jump dy- 
namics between sites and diffusion inside a sphere. 



FLUCTUATING SITE INTERSPACING 

Consider the jump dynamics of a particle between two 
non-equivalent sites separated by a fluctuating distance 
a of stationary correlation function ip(t) and with n and 
T2 being the effective mean residence time of the particle 
in each site. For one site located at the origin, we have 



-PU{r,a) 




, 
a , 



Tl 



otherwise , 



n + t 2 



(10) 



The distribution of relaxation rate is a delta function, 
9n{Qa) = g n [f3s,Qa] = <5„,i with A x = T = rf 1 +t 2 _1 . 
In this case, the normalized scattering intensity is given 
by Eq.(H) with the EISF and the scattering correlation 
function (independent of Q) given by 



A(Q) = (A (Qa)} and C{t) 



-Tt 



(ii) 



where A (Qa) is given in Eq. ([!]). Since T is independent 
of the site interspacing, thus the effect of fluctuations on 
C(t) is simply multiplicative as in dll|). As a result of 
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site interspacing fluctuations, the structure factor A(Q) 
(shown in Fig. [t]) is the EISF for the jump dynamics be- 
tween two sites separated by a fixed distance a averaged 
over the distribution of the distance. As j3e gets smaller, 
the extrema of A(Q) become less pronounced and are 
shifted to the left compare to Aq(Q). The mean-squared 
displacement of the particle is 



where the amplitude 

Ai(Qa) = N- 1 J2n=ijo[2Qasm(n7r/N)]cos(2lmT/N) 
for I = 0, 1, • • • , N, and r is the effective mean residence 
time in each site. Interestingly, the mean-squared dis- 



(X 2 ) = -3 dln[A{Q)] 



d(Q 2 ) 



p(l-p) 



3i? 2 



= p(l-p) (a 2 ) 



R 2 



2(3e 1 + ( 7 r/2/3e) 1 /2 e pe erfc (y/0e) 



(12) 



Solid line in Fig. gB represents this function. When 
j3s — ► oo we have (X 2 ) ~ p(l — p)R 2 like for non- 
fluctuating interspacing while (X 2 ) ~ 3p(l — p)kBT/2<7 
as (3e — > like for motions in a harmonic potential. This 
behavior is compatible and originates from jump dynam- 
ics between soft walls, the broadening of the point sites 
with temperature resulting from the environment fluctu- 
ations. 



A(Q) 




0.3 
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QR 



FIG. 1. The EISF A(Q) in Eq.(|l]J) as a function of QR 
for p = 0.5 and e/keT = 0.5, 2 as quoted on the figure. The 
dashed line represents Ao(Q) in Eq.(^) for a — R and p — 0.5 
for comparison. 



Generalization of Eqs.(|l l|) and ([12|) to the case of jump 
dynamics among TV equivalent sites on a fluctuating circle 
of radius a is straightforward. In this case, the normal- 
ized scattering intensity still writes like in Eq. (^) with the 
EISF given by A(Q) = (A (Qa)) and the Q-dependent 
scattering correlation function by 



C(Q,t) = 



N 



V (MQa)) 



x exp 



2t 



m , (i3) 



placement of the particle is simply 
(a 2 ) is obtained from Eq. (|l2j) . 



'X 1 



i 2 ) where 



FLUCTUATING SPHERE 

The potential for diffusion inside a sphere of radius a 
with impermeable surface is 



U(r,a) 





3C 



r < a , 
r > a . 



(14) 



In this case, the normalized scattering intensity is given 
by Eq.(^) with the scattering correlation function given 
by 



C{Q,t) = {J2 da 



n,l=0 



[A l n (Qa) - A(Q)S i 

[i-^(Q)] 



x e 



2Dt/a2 }m, 

where the A l n (Qa) are given by PJl^l 



(15) 



l \2 



A l n (Qa) 



6(2l + l)(x l n ) 
(xl) 2 - 1(1 + 1) 



Qaji+i(Qa) - lji(Qa) 



(Qaf - (x l n ) 2 



(16) 



~g l n (Qa, a) is the Q-dependcnt effective distribution of re- 
laxation times T l n (a) = a 2 /(x l n ) 2 D, D is the effective dif- 
fusion constant and x l n are the roots of equation p|,p"2||, 
x l n ji + i(x l n ) — lji(x l n ). For the potential in Eq.(||), the 
structure factors (Aq(Q)) (i.e., the EISF for a sphere of 
radius a averaged over the distribution of radii) and the 
EISF of the system A(Q) are therefore: 



(A°(Q)) = (\F(Q)\ 2 ) 
A(Q) = \(F(Q))\ 2 



9^° 



' x 3 [j 1 {xQR)] 2 e-' 3£X dx 



{QR) 2 /~ x 5 e-^ 2 dx 
Sj^x^jtixQR) e'^'dx 



, (17a) 



QR x 5 e~P ex2 dx 



(17b) 



Figure | shows that A(Q) < (A%(Q)) for all QR, and 
the difference between the two functions increases with 
increasing the temperature, i.e. as (3e gets smaller. The 
location of minima of A(Q) and (Aq(Q)) do not coin- 
cide and their values are equal to zero for A(Q) while 
they are different from zero for (Aq(Q)). Like in Fig. [j], 
the location of minima are shifted to the left compare 
to A (Q) and tend to disappear as (3e gets smaller. To 
assess to which extent A(Q) and (Aq(Q)) contribute 
to the elastic intensity, we turn back to Eq.(|l^). The 
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shortest relaxation time of this expansion is T -1 , where 
r = (x\) 2 D/ R 2 is the typical time for a particle to diffuse 
over the entire sphere of radius R. When t ~ T~ l , the 
second term containing the exponential time-dependence 
in Eq.(p^|) can be neglected. In this case the incoherent 
intermediate scattering function reduces to: 

I(Q, if,) = ^(T) (A° (Q)) + [1 - V(T)] A(Q) . (18) 

This expression, which would represent the elastic part 
of S(Q,u>) for the resolution time of about T -1 , con- 
trasts with the case in the absence of fluctuations where 
the incoherent intermediate scattering function relaxes 
to the EISF for time scales of order or greater than 
Figure ||A displays I(Q, tp) versus QR for various values 
of ip(F). Manifestly, I{Q,ip) is different from the EISF 
A(Q) even for ip(T) = 0.1. It is obvious that I{Q,ip) 
will eventually be equal to the EISF in the ^>(r) — > 
limit when the particle diffusion is very slow compare to 
fluctuations of the sphere radius. In the opposite limit 
when the particle diffusion is fast compare to sphere ra- 
dius fluctuations, i.e. ip(T) ~ 1, the I(Q,ip) is essentially 
given by (Aq(Q)). The mean-squared displacement of 
the particle is 



(X 2 ) = -3 



din [7(Q,V0] 



d(Q 2 ) 




2 

1 + 



Q=0 

2 



(3s (f3e) 2 



(19) 



Note that (X 2 ) is independent of ip(T) since, as is illus- 
trated in Fig. 0, in the Gaussian scattering approximation 
the structure factors are, A(Q) = A%(Q) ~ l-3<2 2 (r 2 )/5 
as Q — > 0. Figure |||B shows the temperature dependence 
of the particle mean-squared displacement. At low tem- 
perature (X 2 ) = 3i? 2 /5 like for diffusion inside a sphere 
of radius R while (X 2 ) ~ 9kBT/5cr at higher temperature 
like for motions in a harmonic potential. This behavior is 
compatible and originates from diffusion inside a sphere 
with a thick soft surface, the thickening of the spherical 
surface with temperature resulting from the environment 
fluctuations. 

Let us turn now to the quasielastic term of S(Q,oj). 
The dynamics involved can be characterized by consider- 
ing the position correlation function C (0, t) obtained in 
taking the Q — > limit of C(Q,t) in Eq.(|l5|) to give: 



C(0,t) = 
g((3s,x) 



g((3s, x) e~ rt/x2 dx 



(20) 



x 7 e- 0EX 



J^x 7 e~^ £a;2 dx 



The effective distribution g((3s, x) of relaxation times 
x 2 /V is maximum at x m = max[l, (7/2/3E) 1 / 2 ]. For the 
purpose of illustration we consider the situation where 
the fluctuation dynamics of the sphere radius is described 
by the Langevin oscillator with the correlation function 



= e-^ 2 



cosh(^tt) 



2(i 



sinh(/ii) 



\/7 2 - M 



(21) 



where 7 and u>q are the collision frequency (related to the 
dissipation) and the frequency of the oscillator. The ini- 
tial decay rate constant oftp(t) is zero and its relaxation 
time is j/cOq- 

As shown in Fig. ^A, because of the distribution of 
relaxation times due to fluctuations of the sphere radius 
the C(0, f) is no longer a single exponential. When ip(t) = 
1, the short time behavior of C(0, t) can be fitted by e~ Kt , 
where k is the initial decay rate (see below) of C(0, t), 
while e~ kt , where k is the relaxation rate (see below) of 
C(0, t), is a poor approximation of C(0, t) although they 
coincide around Tt — 15. In the presence of fluctuations 
when ip(t) 7^ 1, C(0, t) is essentially dominated by ip(t)- 

The initial decay rate of C(0, t) is, 



dC(0, t) 



dt 

r 

[3eT/3 



3i? 2 

t =0 5<r 2 ) ' 
(3s -> 00 (low T) 
(3s (high T) 



(22) 



In this example, k is independent of ip(t) and T/k coin- 
cides with the mean-squared displacement. If one uses 
for the diffusion coefficient the relation D = k^T/m^, 
where m is the particle mass and £ the friction coeffi- 
cient, we find that At oc T at low T while k ~ (xi) 2 a/m!; 
independent of temperature for high T. 

The relaxation rate, A: -1 = L C(0,t)dt, which de- 
pends on characteristics of ip(t), is 



r 

k 



dx 



: 136x2 dx 

x 9 (l + -fx 2 /T)e-P £x2 
\+ 1 x 2 /T + {u /T) 2 x A 



(23) 



The ratio k/n as a function of inverse temperature is 
plotted in Fig. [|B. Three types of decays of C(0,t) can 
be noted: k/n < 1 the C(0, t) has a nonexponential de- 
cay dominated by the distribution of relaxation times, 
k/n = 1 then C(0,t) — e~ Kt , and k/n > 1 where the 
nonexponential C(0,t) (which originates from both the 
distribution of relaxation times and il>(t)) may show oscil- 
lations and have long tail decay. These differences in the 
decay of C(0,t) are more pronounced for (3s — > (high 
T) where k/n can be very large while k/n — > 1 (i.e., 
C(0,t) is almost a single exponential) when (3s — * 00 
(low T) . Finally, it goes without saying that when we are 
concern with the jump diffusion in a fluctuating sphere, 
all the expressions derived above remain unchanged ex- 
cept that the exp{— {x l n ) 2 Dt/a 2 } in Eq. (p~5|) is replaced 



byexp j- 



1 _ e -(<b) 2 /2a 2 



jump length and frequency 



jjt\, where b and 7^ are the 
O, respectively. 
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SUMMARY 

Let us summarize the main results derived above. We 
find that as a result of the environment fluctuations the 
EISF of the system is given by Eq.(R) and the incoherent 
scattering correlation function C(t), given in Eq.(j|), is 
C(Q,t) = X; n (e" A " (a)i ) ip(t), where the average is taken 
over the Q-dependent effective distribution of relaxation 
times and rp(t) is the correlation function of the fluctu- 
ating length. 

For the jump dynamics between sites, Eq. (M) reduces 
to (A Q (Qa)) which is the EISF A (Qa) (in Eqtffy) in the 
absence of fluctuations averaged over the distribution of 
site interspacing. As the relaxation rate T (in the absence 
of fluctuations) is independent of the site interspacing, 
the incoherent scattering correlation function C (t) , given 
in Eq. (|Tl|) , is equal to C(t) in the absence of fluctuations 
times ip(t). 

For the diffusion inside a fluctuating sphere, on the 
other hand, Eq.(@) reduces to |(F(Qa))| 2 where F(Qa), 
given in Eq. (J2]) , is the scattered amplitude for diffusion 
inside a sphere of fixed radius a and the average is taken 
over the distribution of radii. However, the normalized 
elastic intensity I(Q,ip), given by Eq.(|18|), involves two 
contributions depending on the relaxation of the sphere 
radius over time scales When the particle diffusion 
is slow compare to the fluctuations of the sphere, I(Q, tp) 
is equal to \(F(Qa))\ 2 while, in the opposite limit, it is 



is given by (\F(Qc 



Since the relaxation time T n {a) 



for the diffusion in a sphere depend on the sphere radius, 
the fluctuations of the sphere radius generate an effective 
distribution of relaxation time so that C{Q,i), in Eq. (15), 

is C(Q,t) = J2n.i=o( e ~ t/T " {a) } where the average 

is taken over the Q-dependent effective distribution of 
relaxation times and if>(t) is the correlation function of 
the sphere radius. 

It worthwhile to mention that when ip(t) = 1, the sit- 
uation resembles to what is known as the heterogeneous 
scenario for the explanation of the stretched exponential 
or Kohlrausch- Williams- Watts relaxation [13]. Such a 
situation may be encountered in the case where the par- 
ticle atoms in a protein system undergo jump dynamics 
between sites (or diffusion inside spherical cages) each of 
them with different site interspacing (radius), i.e. poly- 
dispersity in length scales. In this case, C{Q,t) is expo- 
nential for jump dynamics between sites (provided that 
r is the same for each particle) while it is nonexponen- 
tial for diffusion inside a sphere. The above results with 
tp(t) 7^ 1 are derived in the homogeneous scenario Jl3j 
which assumes that all particles are dynamically identi- 
cal. 

To conclude, we emphasize that the dynamics of par- 
ticle atoms in a protein system is a multidimensional 
problem in which the particles and constituent elements 
of their environment undergo their own dynamics in re- 
spective potentials. The speculative discussion outlined 
above is a quite simplified two-dimensional version (par- 



ticle position r and site interspacing or sphere radius a) 
of the problem. Nonetheless, it is encouraging to see that 
the present analysis indicate that the form of the scat- 
tering intensity, for instance, is quite influenced by the 
fluctuations. This suggests at least to reconsider the way 
of fitting experimental data and revise the interpreta- 
tion of parameters as determined. In the same spirit, an 
analysis combining particle dynamics and fluctuations in 
both the size and shape of the confining geometries can 
be developed along these lines. 
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FIG. 2. The amplitudes (A%(Q)) and A(Q) in Eqs.(p7a|) and ( p7b| ), respectively, as a function of QR for two values of the 
reduced energy e/keT = 0.5, 2 as quoted on the figures. The dashed lines represent Aq(Q) in Eq.(^) for a — R for comparison. 
Note that QR G [0.1, 6] in the typical neutron scattering experiment. 




FIG. 3. Panel A: The intensity I(Q, xji) in Eq.(|is|) for diffusion inside a fluctuating sphere, versus QR for V'(r) = 0.1, 0.5, 0.9 
as quoted on the figure. Panel B: Reduced mean-squared displacements (X 2 )/p(l — p)R? in Eq.(p^) (solid line) for jump 
dynamics and 5(X 2 } /3R 2 in Eq.(|19|) (dashed line) for diffusion, as a function of the reduced temperature ksT/e. 



G 



C(0,t) 




k/K 




l 1 — i — i — n-r 



0.1 



¥ = 1 



_l I I I L_1_U 



1.0 
£/k„T 



_l I I I I I L 



10.0 



FIG. 4. Panel A: Correlation function (7(0,*) in Eq.@ as a function of the reduced time Ft for the reduced energy 
/3e = 0.5. The dashed and dot-dashed lines represent e~ Kt and e~ kt , respectively, with k = 13r/79 and k = 79r/633. Panel 
B: Reduced relaxation rate k/n (ratio of Eqs.(p3|) and (p^)) as a function of e/keT. For the two panels, we have wo = T/2 
and the quotations 'under' (7 = r/4V 'crit'(7 = V), and 'over' (7 = 4F) correspond respectively to underdamped, critical and 
overdamped regimes of ip(t) in Eq.(pil) 
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